On the black hole singularity issue in loop quantum gravity 
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Abstract 

This paper presents a brief overview on the issue of singularity resolution in loop quantum gravity presented at the 
Theory Canada IV conference at the Centre de Recherches Mathematiques at the Universite de Montreal (June 4-7, 
2008). The intended audience is theoretical physicists who are non-specialist in the field and therefore much of the 
technical detail is omitted here. Instead, a brief review of loop quantum gravity is presented, followed by a survey 
of previous and current work on results concerning the resolution of the classical black hole singularity within loop 
quantum gravity. 

Resume 

Nous presentons une breve revue portant sur la resolution du probleme de singularite en gravite quantique a boucles, 
destinee a la conference Theorie Canada 4 tenue au Centre de Recherches Mathematiques de l'Universite de Montreal 
(4-7 juin 2008). L'audience visee est celle de physiciens theoriques qui ne sont pas specialistes du domaine et par 
consequent, nous omettons la plupart des details techniques. A la place, nous prsentons une breve revue de la gravite 
quantique a boucles, suivie dun survol des resultats anterieurs et presents concernant la resolution de la singularite 
classique du trou noir dans le cadre de la gravite quantique a boucles. 
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1 Introduction 

General relativity is arguably the most successful theory 
of gravity to date. Even after more than 90 years, the 
theory is extremely robust at predicting solar system scale 
physics to astounding accuracy. On much larger scales, 
general relativity supplemented by various exotic sources, 
such as dark matter or dark energy seems to explain very 
well the properties of the universe we live in. However, on 
the extremely small scale, or equivalently very high energy 
scales, there is precious little we know about gravity. In 
this sector we have to heavily rely on what we expect grav- 
ity to behave like in this regime. At very high energies 
we expect that classical general relativity will get modified 
in some way by quantum effects. However, the quantiza- 
tion of the gravitational field is not at all straight-forward, 
as full diffeomorphism invariance requires a different ap- 
proach than those applied to conventional field theories. 
Other problems associated with the quantization of grav- 
ity may be found in pQ, which also presents a nice survey 
of the field in general. 



One potential candidate for a theory of quantum grav- 
ity is known as loop quantum gravity which, over the past 
couple of decades, has become a strong contender in this 
arena. The loop quantum gravity program incorporates 
diffeomorphism invariance as demanded by general rela- 
tivity and, as will be outlined below, is in fact a quanti- 
zation of general relativity. An excellent (but technical) 
overview of loop quantum gravity may be found in [2]. 
Essentially, what arises out of loop quantum gravity is a 
theory of quantum geometry, which is quite appealing from 
the point of view of general relativity. 

Having a candidate theory at hand, the natural question 
to ask is what theoretical (in the absence of experimental) 
issues do we wish the theory to resolve? There are proba- 
bly many but two that come to mind are the explanation 
of black hole entropy and the elimination of the classical 
singularities such as those encountered in cosmology and 
black hole physics. 

Regarding the entropy, there is still debate on the actual 
source of this entropy. One belief is that the source is 
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strictly gravitational in origin. That is, one should be able 
to define microstates in a full quantum theory of gravity 
which, when counted, should yield the correct entropy law. 
It has been shown that for spherical black holes the correct 
entropy law is derived within loop quantum gravity [3] - 
[14] . The result has been extended to higher genus black 
holes, where it was found that the sub-leading term to 
the entropy is genus dependent [15 . This is in agreement 
with studies of black hole entropy of higher genus black 
holes that have utilised non-quantum gravity techniques 
[16], [17] [18]. Progress from the spin-foam approach has 
also been made [19], [20] , 

On the singularity issue, most of the work in loop quan- 
tum gravity has been done within the paradigm of loop 
quantum cosmology (LQC). In LQC all but the classi- 
cally relevant degrees of freedom are frozen out first and 
then the resulting, symmetry reduced, system is quantized. 
There are, of course, some limitations inherent in the var- 
ious approaches to LQC (for example, see [21]). However, 
studies to date indicate that the quantum corrections be- 
come important as one approaches the big bang singularity. 
The singularity itself is therefore replaced with a smooth 
bounce. 

There are many studies in the literature on loop quan- 
tum cosmology. We do not survey the vast literature here 
due to space limitations. However, the interested reader is 
referred to [22] and the references therein. Instead, here 
we chose to briefly discuss the progress on the black hole 
front, which utilises many of the same techniques used in 
LQC. 

It should be stressed that the intended audience is not 
the expert, the Theory Canada conferences being made up 
of theoretical physicists from a variety of fields and back- 
grounds. Instead, only the main points are presented, with 
little technical detail. However, it is hoped that this will 
give the reader an overview of loop quantum gravity as well 
as the flavor of the various approaches to the problem of 
singularities within this theory. Of current research inter- 
est to us is the extension of these studies to various other 
singular solutions of general relativity. Progress on this 
front has been made in |23] where singularities have been 
shown to be avoided in black holes other than spherical. 



where g is the determinant of the metric tensor, g^ v , and 
i?(g M „) is the Ricci scalar, which is to be viewed as a func- 
tional of the configuration variable, which is the metric 
tensor. The quantity Im{9iiv) represents the matter field 
action. Demanding that the variation of this action with 
respect to g^ vanishes gives rise to the Einstein field equa- 
tions, 

Rpu — -^Rg^v — 87rT M ^ , (2) 

with R^ v being the Ricci tensor and T^ v the stress-energy 
tensor of the matter field. 

We will consider here the canonical quantization of the 
Einstein-Hilbert action within the paradigm of loop quan- 
tum gravity. A covariant form also exists, giving rise to 
spin-foam models, so named due to the resemblance of the 
Feynman diagram analogs to a foam of bubbles. However, 
for an introduction the canonical approach is much more 
appropriate. In the canonical approach one splits space- 
time into "space" and "time" by writing the line element 
associated with the metric as: 



ds 2 



g^dx^dx" 



(N 2 - N a N a ) dt 2 
2N b dx h dt + q ab dx a dx b . 



(3) 



Here, N and N b are the lapse function and shift vectors 
respectively and q ab is the metric on the three-dimensional 
hyper-surfaces, S. These quantities are shown schemati- 
cally in figure [T] 




ds 3 =-(N -N tf Wf 2 + 2Ndx b dt +q dxdx 



Figure 1: ADM decomposition of space-time into "space" and 
"time". In this decomposition, one studies how the three ge- 
ometry of £ evolves with time t. 

In the variables of equation ([3]), the action (fTJ) may be 
written as 

1 



2 A brief introduction to loop quantum gravity j _ 1 i dt I dx 3 Tl ab a +2A^V^ ( a~ 1 ^ 2 U ab 

We begin by considering the Einstein-Hilbert action of gen- 
eral relativity: 



+Nq 1 ' 2 [ - q-'n^ + \q^ 2 



(4) 



1 

8tt 



d 4 x y^gR(g^) + hi(g^) ■ 



with n := U ab q ab , n ab := q- 1 ' 2 [K ab - K q ab ] , K := 
(1) K ab q ab . K ab is the extrinsic curvature tensor, the 
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Ricci scalar of the three-dimensional hyper-surface and q 
the determinant of q ab with over-dots denoting differentia- 
tion with respect to the slicing time, t. In this scheme, the 
field variable is q ab and its conjugate momentum is IP 6 . 
The resulting equations of motion are: 



-l/2j-ph 



3 



1/2 



RW-q-'n^ + lq-W 



= 0=: ~V b , 
= 0=:-3. 



The action in (HJ) leads to a Hamiltonian density: 
n G = -^-U ab q ab -£ G = N b V b + NS, 

with the following symplectic structure: 

{n ab (x), q cd (y)} = 16n61 c 6%S( X , y), 
|n Qb (x), n cd (y)} = = {q ab (x), q cd (y)} . 

After quantization, this leads to the well-known Wheeler- 
DeWitt equation: 



1 

— qabqcd + -^lacQbd 



_S 6_ 

Sq ac Sqbd 



+ qR i3) \v(q) = 0. (5) 



The above equation was first derived by DeWitt [24j and 
Wheeler [25] in the 1960s. There are some problems as- 
sociated with the approach leading to (|5|). Some of the 
problems are associated with all canonical approaches to 
quantum gravity, such as the "issue of time" , which leaves 
open exactly how to determine evolutions in a theory where 
time is not manifest in the quantum equations. Other is- 
sues with the Wheeler-deWitt formalism involve inconsis- 
tencies in some transition amplitudes in a path integral 
formulation of the theory. Also, the three-metric config- 
uration variable does not appear as a gauge field in the 
theory, which means many of the techniques used in gauge 
field theories are not applicable here. As well, as will be 
briefly mentioned below, the singularities present in clas- 
sical gravitation are not eliminated within the WdW for- 
malism. 

In order to make contact with gauge theories, Ashtekar 
formulated the gravitational action in terms of densitized 
triads E°i [26], which are related to the 3- metric compo- 
nents via 

qq ab = EIE^ , 
and are related to the standard triad, e\, by 



E a — -e abc e- iV p k 



In terms of these new variables, the 3+1 Einstein-Hubert 
action may be cast in the form: 



I= h h " 



E a 4 K\ ~N h V b - NS 



-N i e ijk E a3 K k a ] d 3 



K\ := 



v/dctp) J 



The canonically conjugate variables are now E® and K l a . 
The transformation group of the triad in E is 5*0(3) and 
as such, so is that of its conjugate momentum, K l a . The 
natural 50(3) connection compatible with the vanishing 
of the covariant derivative of the triad is given by T l a as in 



\aC\] + 6 jfc r [a e b] 



= 0. 



(6) 



We can also define another so(3) connection A l a by the 
linear combination: 

A* a : =1 K\ l+ r a , 

where 7 is known as the Immirzi parameter, which must 
be determined in some way and is therefore a type of 
quantization ambiguity at this point. The new connec- 
tion A % a turns out to also be conjugate to E % a and this 
pair, (A l a , E 1 ^, make up the Ashtekar-Barbero variables 
on which loop quantum gravity is based. In terms of these 
new variables, the action is written as 



'-h 



dt 



E a ;A\ 



N b V b - 



NS - N l G. 



d 6 x 



Gi :=d a E a l + e l3 k Ai a E a k 
Here. 



H:=^ b -(l+ 7 2 )if i b G i , 



(7) 



(8a) 



5:=[det(£)]~ 1/2 E 1 E 



fc k r ab 



2(l + 7 2 )K\ a K{ 



(8b) 



with F\ 



d a A\ ~ d b A\ 



i A3 A k 
jk A a A b ■ 



The corresponding symplectic structure is given by: 

{25».(x), A* 6 (y)} =8TTjd a b 6 i j 6(x,y), 
{^(x), E) (y)} = {A* a (x), A\(y)} 



0. 



To the above system one could also add matter couplings 
by supplementing the action with a matter term and quan- 
tizing appropriately. The problems associated with ([5]) are 
not present in this representation and it can be seen that 
A % b enters like a gauge field. 
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Variation of the action with respect to the shift-vector, 
iVft, and lapse function, N, supplemented with the fixing of 
rotational freedom leads to what are known as the vector, 
scalar and Gauss constraints respectively: 



V b =0, 

s=o, 

Gi =0 



(9a) 
(9b) 
(9c) 



(the explicit forms are given below.) Note that now there is 
a constraint equation for Gi (the Gauss constraint). This 
reflects the fact the triad possesses a rotational freedom; 
one can choose different frames locally by rotating the 
triad. This redundancy is eliminated in the new Gauss 
constraint. 

It is customary to work in SU(2), the universal covering 
group, since the Lie algebra is the same and this change 
does not affect the constraints. It is also useful from the 
point of view of bridging quantum gravity with standard 
gauge field theory. Therefore the connection A is now to 
be viewed as an SU (2) connection. 

The constraint equations that arise, plus the rotational 
degree of freedom of the triad give rise to the following 
equations after quantization: 



Gi |*) =A^ |*) = 0, 



sm = 



-.E°iE b , 



(10a) 
(10b) 



det(£) 



e ij k F k ab - 2(1 + 7 2 )K\ a K j b] ] |*> = , (10c) 

where D a is the gauge covariant derivative given in the 
second line of equation {?])• 

One can construct functionals ^(A) which are function- 
al of A and are square-integrable with respect to the 
Ashtekar-Lewandowski-Baez measure 0: 



dfi[A] V[A] $L4] 



This measure is gauge as well as diffeomorphism invariant. 
The corresponding Hilbert space possesses some peculiar 
properties which we do not address here. The interested 
reader may refer to [2] and references therein. 

Before continuing, we summarize a few points about loop 
quantum gravity that the reader may be curious about 



1 Strictly speaking it is the holonomies of A which are of funda- 
mental importance. However, as we are omitting technical details, 
this can be ignored at this stage. It is important to note however that 
with the concept of holonomy, paths become an important feature of 
the theory. 



(recalling that the intended audience is a theorist who does 
not specialize in loop quantum gravity): 

1. The scheme is background independent and respects 
diffeomorphism invariance. The choice of time slicing 
is arbitrary and does not affect the physics. 

2. A superpartner can be accommodated and therefore 
supersymmetry can be incorporated. This has been 
done [27] . 

3. Instead of the Einstein-Hilbert action, one can accom- 
modate geometric actions made up of arbitrary cur- 
vature invariants. The scheme is generally similar to 
that outlined above. 

4. Higher dimensions can be accommodated. 

5. The theory is divergence free. 

It should be noted that 2, 3, and 4 are not required but 
simply can be accommodated. 

2.1 Geometric operators 

Having discussed (very briefly) the loop quantum gravity 
program we now turn our attention to important operators 
within the theory. As this is a theory of quantum geometry, 
areas and volumes are of particular importance. "Classi- 
cal" areas (Ar) and volumes (V) may be constructed from 
the triad and, equivalently, from the densitized triad: 



Ar = J dfa \/E a i EJn a n b , 



V 



f d 3 x [\e abB f* k E a i E b j E%\] 
Jm 



Since these involve E, they become operators in the quan- 
tum theory. The area operator possesses the following 
eigenvalues: 



(Ar) |¥> = 8nl 2 p7 ^j(j + 1) |*) 



(11) 



l p being the Planck length. The interpretation of the above 
equation is as follows; The state functions of the theory 
describe a spin-network. This arises from the fact that 
the relevant operator made up of the configuration vari- 
able is not the connection A itself, but the holonomy of A 
along an oriented path e. This holonomy is a functional of 
the connection A and allows for the parallel transport of 
spinors along e. Combinations of these paths make up the 
spin-network and each edge (e) of the path carries with it a 
half-integer spin, denoted by the symbol j in (fTT|) . When a 
segment of the spin network penetrates a surface at a point 
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p, it endows the surface with an amount of area related to 
its spin. The specific amount is given by the eigenvalue in 
(fTTj) . The total area of the surface is given by adding all of 
the eigenvalues of the penetrating network segments over 
the surface S. A simple diagram illustrating this is shown 
in figure [50 




Figure 2: A spin-network "puncturing" a surface. Each part of 
the spin-network puncturing the surface endows it with a cer- 
tain area and local angular deficit, yielding non-trivial geometry 
for the surface. The nodes (n) are associated with volumes. 



The very interesting issue is that this theory predicts 
that, on the small scale, space is discrete\ ft should be 
emphasized that this is a prediction of the theory and is 
not put in "by hand" . The volume operator is more com- 
plicated to deal with but also predicts that volumes are 
quantized. (Volumes are associated with the nodes of the 
spin- network.) 



first imposes certain symmetries at the classical level and 
then quantizes the symmetry reduced model in the hope 
that the essential features of the full theory are retained. 

There are several approaches to studying black hole sin- 
gularities in loop quantum gravity. We will first be con- 
centrating on the approach used in |28j (related works may 
be found in [5^ and [3D] and references therein.) 

To date, the Schwarzschild black hole has received by 
far the most attention. The first step is the symmetry 
reduction to a symmetric system appropriate for the black 
hole interior (Rx S 2 ). The connection and triad respecting 
this symmetry are given by 

A =A r t 3 dR + (Am + Am) d9 + [(Am - Am) sin9 
+r 3 cos 6] dip , (12a) 



E =E r t 3 sin( 



+ (EWi +E 2 t 2 ) sin (9 



(E 1 



T2 



dR 



00 



d_ 

dip ' 



(12b) 



where the Tj denote the standard su(2) basis. (Note that 
the indices 1 and 2 on the quantities A and E are not di- 
rectly related to the 1 and 2 indices on the r matrices.) 
In this reduced scenario, the vector constraint is auto- 
matically satisfied and the Gauss contraint is satisfied if 
A\E 2 — A 2 E X = 0. In the sequel A\ and E 1 are set equal 
to zero to satisfy this constraint. The Hamiltonian con- 
straint, which is a combination of the Gauss, vector and 
scalar constraints yields: 



H = - 



N sgnE R 
&H 2 y^E^E 2 



[{A 2 2+1 2 ){E 2 ) 2 + 2A R E R A 2 E 2 } . 

(13) 

It can be shown (still dealing at the classical level), that 
the following quantities solve the equations of motion [28 : 



O Singularity resolution 

It has long been thought that any viable theory of quantum 
gravity should eliminate the singularities predicted by the 
classical theory. Specifically, the space-like singularities of 
black hole physics and standard cosmology should be re- 
moved. If a theory is successful in this endeavor, a natural 
question to ask is what replaces the singularity? We will 
briefly outline progress and results within loop quantum 
gravity on this front here. 

Given the complexity of the situation, studies are gener- 
ally carried out in symmetry reduced models. That is, one 



2 We are making an assumption here regarding how the spin- 
network pierces the surface S. The general case yields eigenvalues 
which arc slightly more complicated than (lilt . 



A 2 = ± 7v /(2Af-T)/T, E 2 — Eq\/T(2M - T) 
Ar= TA° r1 ^, E r = ±T 2 . 

(A coordinate re-scaling has tacitly been assumed.) 

It turns out (not surprisingly) that the above (densi- 
tized) triad is compatible with this line element's metric: 

d ° 2 = ~ A d R 2 + T 2 d0 ' 2 + T 2 si n 2 dip 2 , 

which is the T-domain, or black hole interior, of the 
Schwarzschild space-time. 

Next, the symmetry reduced system is quantized via: 



E 2 := -i^l 



d 
dAo 



E R := -21^1 







dA R ' 
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with the following spectrum: 



be shown that 



E 2 |M, t) = plfr \n, t) , E R r) = 7 Z*r | M , T ) . 

Here, t and /1 represent the coordinate length of the holon- 
omy curve in the R and angular directions respectively in 
the state vectors. 

The volume operator is also constructed: 



V = 4ir\E 2 \^\E R \ , 
with eigenvalues v = 27r-y 3 / 2 Zp v \r\ 



(14a) 
(14b) 



Also, in what follows, a general state will be expanded in 
eigenfunctions of the triad as: 



so that it is sufficient to study effects on t) and on the 
expansion coefficients C^t- 

Before proceeding we should note that there are several 
ways to indicate potential singularity avoidance. One is to 
study the quantum analogues of the curvature. Another is 
to study the evolution equations and see if they are smooth 
(non-singular) where the classical singularity exists (T = 
0). Both were analyzed in [25], whose results we quote 
here. For the curvature argument, we require the co-triad, 
u: 

lu = LJRT3 dR + W2T2 d9 — UJ2T1 sin 6 dip . 

Classically, the co-triad involves inverse powers of the triad 
and is useful for forming quantities measuring curvature 
and in the Hamiltonian constraint, which also contains in- 
verse powers of the triad. Explicitly, the classical co-triad 
is given by: 



ljr = sgn (E R ) 



\E 2 



lo 2 = sgn (E 2 ) J \E R \. (15) 



To calculate the quantum quantities appearing in the cor- 
responding operator, ujr, which is the component involving 
the inverse power of the triad, one first rewrites the classi- 
cal co-triad in a classically equivalent manner and creates 
the corresponding quantum operator by analogy [31: 



1 



27r 7 /2 



Tr{T 3 h R [h-\ F]} 



(16) 



Here Hr is the holonomy (operator) of Ar, which has a 
well defined inverse operator. With this definition it can 



4^7/2 
V7, 



(v^t+i ~ v^t-i) |m, t) (17) 



= ^im(V¥^\-V^i\)\»,t) , 



UR \H, T) 



&2 |Mi t ) = v^pSgnOw l r l Im, r) . 

Note that the C02 eigenvalues have the same relation to 
E s eigenvalues as in the classical relationship. However, 
u) R 's eigenvalues relation to the triad's eigenvalues deviates 
from classical relation and only asymptotically approaches 
the classical relation for large r. Since the classical singu- 
larity is located at t = 0, quantum effects become impor- 
tant near the singularity, ft may be seen from the above 
that the co-triad is not singular in the quantum theory and 
t = is well behaved. 

The evolution of the Hamiltonian constraint has also 
been used to argue the elimination of the classical singular- 
ity |28j . Utilising the quantum operators, the Hamiltonian 
constraint can be constructed. The Hamiltonian constraint 
operator, acting on the states yields a difference equation 
for C^r of the form [28] : 

tff tav C,. r = 26 (V|r + 2<J| + vV\) [C„ + 2s, T + 2s 
-C„ 25, t -\- 2s} + (V\r + 5\ - VF~6\) [fa + 2S)C„ + 4 ,, r 



-(1 + 2 1 2 S 2 ) f iC,, T + (n - 25)C^u,r) 



+ 28(^\t~25\ + ^\t~\) 



■ 25, t — 25 + 26, t - 26 



= 0. 
(18) 



The quantity 5 is sometimes known as the holonomy pa- 
rameter, which measures the size of the steps taken in the 
R and angular directions (5 is assumed to be the same in 
all directions here.) 

The particular expression (|f 8p is valid for /i > 45 but 
similar expressions exist for other values of /x. The form of 
the constraint for /1 < 4 dictates that only data at [i = 
and /i = S needs to be specified and the issue of singularity 
resolution turns out to be insensitive to the choice of initial 
data [28]. 

The evolution described by ([T8]) procedes smoothly as 
long as the coefficient of the wavefunction is non-zero at 
t — 25. The coefficient here is y/\r — 25\ + y^jrf and never 
vanishes. The evolution is therefore smooth across the sin- 
gularity. It turns out that the analogous Wheeler-deWitt 
equation cannot be evolved beyond the point correspond- 
ing to the classical singularity. 

Before proceeding it is worth noting that some issues ex- 
ist which make the above analysis potentially suspect. It 
has been argued that by considering 6 to be a constant the 
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above system does not possess an good semi-classical limit 
[32] . In the field of loop quantum cosmology, where similar 
studies are performed, the parameter analogous to S must 
depend explicitly on the scale factor to get an appropri- 
ate semi-classical limit. Therefore, in [33], the condition of 
constant S was relaxed. The approach also relied on an ef- 
fective Hamiltonian approach, whose equations of motion 
capture the essence of the main quantum corrections. This 
is more in-line with approaches taken in loop quantum cos- 
mology and was first applied to black holes in [34] for the 
case of constant 6. For constant 5, the effective dynamics 
are encoded in the modified Hamiltonian (compare with 

my- 



yields 



Keff - - 



N 



2 sin{A 2 5) sin(A R S) ^^ 



sin A 2 S n 



E 2 



(19) 



An appropriate choice of lapse function, N, allows one 
to solve the resulting system of Hamilton equations for 
Ar, A 2 , E r and E 2 . The mathematical result may be 
found in [33] and [33J . Qualitatively, the solutions describe 
a bounce in place of the classical singularity. This bounce 
connects two (non-singular) black holes of generally un- 
equal mass. 

The above results, though promising, depend on a fixed 
holonomy parameter, which has been shown to be problem- 
atic in the semi-classical limit when similar fixed 5 tech- 
niques were used in loop quantum cosmology. Therefore 
an improved scheme was sought in [33] . Following loop 
quantum cosmology, non fixed S quantities are utilized in- 
stead. The natural question that arises then is how to set 
the fundamental size of the two 5 parameters that appear 
in the improved theory. We denote the holonomy param- 
eter associated with the "radial" direction as Sr and the 
holonomy parameter associated with the angular direction 
as 5 2 , to make ties with the notation above. 

The holonomies in the Hamiltonian involve closed loops 
and one may constrain the size of the loop by appealing to 
the fact that the holonomy loops make up the minimum 
area allowed by loop quantum gravity 0: 



Ar = iVSirl; 



7- 



(20) 



The area of a 9 — <f> surface element is given by 

Ar e ^E R (6 2 f , (21) 
which, by restricting that this equal the minimum area, 



3 There are various schemes to impose a non- fixed <5 and we only 
concentrate on one here for brevity. 



Aro 
E R 



(22) 



Similarly, given that the area of an R — 9 surface element 
is given by Atrq = E 2 6rS 2 , one arrives at the condition 



y/Ar E R 
E 2 



(23) 



The generalization of ([19]) to form the Hamiltonian con- 
straint is straight-forward: 



N 



5 2 Sr 



sin 2 A 2 5 2 2 
(S 2 ) 2 +7 



E 2 



(24) 



The resulting equations of motion were solved numerically 
in [33] , There, the authors started the evolution near the 
horizon, using initial data compatible with the classical 
solution, along with the Hamiltonian constraint TL = 0. As 
in the previous scenarios, the singularity is avoided, with 
2-sphere radii being bounded from below. In this improved 
scheme the solution does not connect two black holes but 
instead asymptotes to a Nariai solution at late time. For 
reference, the type of Nariai solution compatible with this 
evolution takes the form 



ds 2 = a 2 (dr 2 + cosh 2 r d X 2 ) + b 2 dfl 2 . 



(25) 



with a and b constants. Several figures are presented below 
in figure [3] to show a sample evolution (see figure caption 
for details). 



4 Concluding remarks 

We have briefly summarized here the status of black hole 
singularity studies within loop quantum gravity. It should 
be noted that at this stage it is not clear if the limita- 
tions of the approximations above are serious enough to 
disqualify all of the results. A discussion of the drawbacks 
and strengths to various methods within the framework of 
loop quantum cosmology may be found in the recent pa- 
per [35]. However, it is encouraging that all the methods 
indicate that the classical singularity is indeed avoided in 
the evolution. One may say with some confidence that the 
singularity may therefore be avoided within the full theory. 

Future studies of interest would be to study other sys- 
tems which are singular in classical gravity as well as the 
incorporation of higher order effects. Of particular interest 
in all of these studies is what exactly replaces the classi- 
cal singularity. This is certainly one of the most engaging 
questions a quantum gravity theory can answer. 
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Figure 3: A sample evolution utilizing the improved quantiza- 
tion technique in [33] . The top plot relates the two densitized 
triad variables, E R and E 2 . At the classical singularity, the 
curves would intersect the origin. This does not happen in the 
quantum evolution. On the bottom, the combination A R /E 2 as 
a function of the Hamiltonian evolution time t (T = e*). This 
remains finite and approaches a constant for t — > — oo. This 
quantity would become unbounded at a curvature singularity. 
I thank Steve Kloster for generating the figures. 



In closing I would like to thank Richard MacKenzie, 
Arundhati Dasgupta, and all the other personnel involved 
with the organization of the Theory Canada IV conference. 
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